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Abstract 

We investigate non-degenerate Lagrangians of the form 

f(u x ,u y ,u t ) dx dydt 

such that the corresponding Euler-Lagrange equations (f Ux )x + (fu v )y + (fu t )t = are 
integrable by the method of hydrodynamic reductions. We demonstrate that the integra- 
bility conditions, which constitute an involutive over-determined system of fourth order 
PDEs for the Lagrangian density /, are invariant under a 20-parameter group of Lie-point 
symmetries whose action on the moduli space of integrable Lagrangians has an open orbit. 
Lhe density of the 'master-Lagrangian' corresponding to this orbit is shown to be a mod- 
ular form in three variables defined on a complex hyperbolic ball. We demonstrate how 
the knowledge of the symmetry group allows one to linearise the integrability conditions. 

MSC: 35Q58, 37K05, 37K10, 37K25. 
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1 Introduction 



In this paper we investigate integrable three-dimensional Euler-Lagrange equations, 



(fu x )x + (fu y )y + (fu t )t 



0. 



corresponding to Lagrangian densities of the form f(u x ,u y ,Ut). Familiar examples include the 



dispersionless KP equation u 



xt 



Ua*Un 



u 



uu 



with the Lagrangian density / 



lul+ul~ //.,///: 



this equation, also known as the Khokhlov-Zabolotskaya equation, arises in non-linear acoustics 
[19|. Another example, u xx + u yy = e Ut u t t, is known as the Boyer-Finley equation pj: it appears 
as a symmetry reduction of the self-duality equations, and corresponds to the Lagrangian 
density / = u\ + u 2 y — 2e Ut . 

The paper [8] provides a system of partial differential equations for the Lagrangian density 



f(a, b, c) (we set a 



u , 



u 



c = u t ) which are necessary and sufficient for the integrability 



of the equation (11.11) by the method of hydrodynamic reductions as proposed in [7J. These 
conditions can be represented in a remarkable compact form: 

Theorem 1 [8]. For a non- degenerate Lagrangian, the Euler-Lagrange equation U.l\) is 
integrable by the method of hydrodynamic reductions if and only if the density f satisfies the 
relation 

d 4 f = d 3 f w + -det(dM); (1.2) 

here d 3 f and d 4 f are the symmetric differentials of /. The Hessian H and the 4x4 matrix M 
are defined as follows: 



H = det 



f aa fab fac 
fab fbb fbc 
fac fbc fee 



M 



( f a f b f c \ 

fa f aa fab fac 

fb fab fbb fbc 

\ fc fac fbc fee ) 



;i.3) 



The differential dM = M a da + M^db + M c dc is a matrix- valued form 



/ 

faa 
fab 
\ fac 



faa 
faaa 
faab 
faac 



fab 
faab 
fabb 



fac ' 
faac 
fabc 



da + 



fabc face ) 



( o f ab 

fab faab 

fbb fabb 

\ fbc fabc 



fbb 
fabb 
fbbb 
fbbc 



fbc \ 

fabc 
fbbc 
fbcc / 



db + 



( 

fac 
fbc 
\ fee 



fac 
faac 
fabc 
face 



fbc fee ^ 

fabc face 

fbbc fbcc 

fbcc f ccc ) 



dc. 







A Lagrangian is said to be non-degenerate iff H ^ (we point out that the equations H 
and detM = have been discussed in the literature, see [6] and references therein). 

Both sides of the relation (ll.2p are homogeneous symmetric quartics in da, db, dc. Equating 
similar terms we obtain expressions for all fourth order partial derivatives of the density / in 
terms of its second and third order derivatives (15 equations altogether). The resulting over- 
determined system for / is in involution, and its solution space is 20-dimensional: indeed, the 
values of partial derivatives of / up to order 3 at a point (ao,&o,co) amount to 20 arbitrary 
constants. Thus, we are dealing with a 20-dimensional moduli space of integrable Lagrangians. 
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In Sect. 2 we prove that the integrability conditions (11.21) are invariant under a 20-parameter 
group of Lie-point symmetries whose action on the moduli space of integrable Lagrangians 
possesses an open orbit. 

Explicit formulae for integrable Lagrangians in terms of modular forms are constructed in 
Sect. 3. We first consider Lagrangian densities of the form / = u x u y g(u t ), which can be viewed 
as a deformation of the integrable density / = u x u y u t found in [8] . By virtue of the integrability 
conditions (ll.2p . the function g has to satisfy the fourth order ODE 

g""(g 2 g" - 2g(gf) - 9(</) V) 2 + 2gg'g"g"' + WW ~ ^'"f = 0, 

which inherits a remarkable Gl(2, M)-invariance. We prove that he 'generic' solution of this 
ODE is given by the series 

g (ut) = e (a2 - a/3+/32)ut = 1 + 6e Ut + 6e 3ut + 6e 4ut + 12e 7ut + ....; 

(a,/3)eZ 2 

notice that under the substitution w t = 2wiz the right hand side of this formula becomes a 
special modular form of weight one and level three, known as the Eisenstein series Ei^(z). 
We point out that modular forms and non-linear ODEs related to them appear in a variety of 
problems in mathematical physics, see e.g. [TJ [2], SI [9j [101 [HI [151 [18] and references therein. 

Lagrangian densities of the form g{u x ,u y )u t and the general case f(u x ,u y ,u t ) are discussed 
in Sect. 3.2 and 3.3, respectively. Here the 'generic' solution is an automorphic form of two 
(three) variables. 



2 Symmetry group of the problem 

The first main observation, overlooked in [8], is the invariance of the integrability conditions 
(11.21) under projective transformations of the form 

_ _ h(a,b,c) ~ _ l 2 (a,b,c) _ _ h(a,b,c) ~ _ f 

l(a,b, c) ' l(a,b, c) ' l(a,b, c) ' l(a, b, c) ' 

here l,h,h,h are arbitrary (inhomogeneous) linear forms in a,b,c. Introducing the quartic 
form 

F = Hd 4 f - d 3 fdH - 3det(dM), 

one can verify that F = l 4 F, which establishes the SL(4, M) -invariance of the integrability 
conditions (II. 2p . Combined with obvious symmetries of the form 

/ = s f + aa + (3b + 7c + 5, (2.5) 

this provides a 20-dimensional symmetry group of the problem. 

Remark. The class of Euler-Lagrange equations (11.11) is form-invariant under a point group 
generated by arbitrary linear transformations of the variables x, y, t and u. Obviously, point 



3 



transformations preserve the integrability. Since the prolongation of these transformations 
to the variables a,b,c and / is given by (12. 4p . this explains the SL(4, M)-invariance of the 
integrability conditions (11. 2p . 

The main result of this section is the following 

Theorem 2. The action of the symmetry group on the 20- dimensional moduli space of 
integrable Lagrangians possesses an open orbit. 

Proof: 

The infinitesimal generators of the symmetry group ( 12.41) . ( 12.51) are the following vector fields: 
3 translations in a, b, c: 



d d d 
da' <% ' dc 



9 linear transformations of a, b, c: 



d d d d d d d d d 

a a~' c «~' b ^ a ^~> c a~> 

oa oa oa ob ob ob oc oc oc 

3 projective transformations of a, b, c, f: 

2 d . d d .8 . d l2 d . d . . d d , d 2 d .8 

a — + ab— + ac— + a/—, a6— + b — + be— + bf — , ac— + be— + c — + cf — ; 
aa ao ac of oa ob oc of oa ob oc of 

moreover, we have 5 extra generators corresponding to the transformations (12.51) : 

d 3 d d d 

df' a Tf df' c d~f } df' 

The main idea of the proof is to prolong these infinitesimal generators to the 20- dimensional 
moduli space of solutions of the involutive system ( 11. 2ft . We point out that, since all fourth 
order derivatives of / are explicitly known, this moduli space can be identified with the values 
of / and its partial derivatives f i} fy, up to order three (20 parameters altogether). The 
prolongation can be calculated as follows: 

(1) Following the standard notation adopted in the symmetry analysis of differential equations 
[TTI [16], we introduce the variables x l = a, x 2 = b, x 3 = c and represent each of the above 
generators in the form 

i_d_ + d_. 
dx i ^ df ' 

here £ l and r] are functions of x l and /. 

(2) Prolong infinitesimal generators to the third order jet space with coordinates x\ f, /j, fy, f^, 

i 9 d d d d 

^ + V df +Q dK + Qj df~j + Qik dh~k 
where Q, Qj and Qjk are calculated according to the standard prolongation formulae 

Ci = Di{n) - fsDiiC), Co" = DjCi - fisDjiC), Qjk = D k Cij - f ija D h (Z')- (2.6) 
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here Di denotes the operator of total differentiation with respect to x l . 

(3) To eliminate the ^-terms we subtract the linear combination of total derivatives t^Di from 
the prolonged operators where, in Di, it is sufficient to keep only the following terms: 





Ox' -"Of ■'■'Of, ■"'<)!,, r7/ ; 



kl 



notice that, since fijki are explicit functions of lower order derivatives, the resulting opera- 
tors will be well-defined vector fields on the 20-dimensional moduli space with coordinates 
fifiifijifijk- Although these operators will depend on the variables x % as on parameters (in- 
deed, the isomorphism of the moduli space with the space /, f, fy, fyi- depends on the choice 
of a point in the x-space), all algebraic properties of these operators will be x- independent. 

(4) Finally, the dimension of the maximal orbit equals the rank of the 20 x 20 matrix of 
coefficients of these operators. It remains to point out that this rank equals 20 for any 'random' 
numerical choice of the values for x\ f, f, fj, 

3 Lagrangian densities in terms of modular forms 

In this section we provide explicit formulae for integrable Lagrangians in terms of modular 
forms. We start with the case f(u x ,u y ,u t ) = u x u y g{u t ) (Sect. 3.1), where g is shown to 
be an Eisenstein series E\^\ a special modular form of weight 1 and level three. The case 
f(u x ,u y ,u t ) = g(u x ,u y )u t and the general case f(u x ,u y ,u t ) are discussed in Sect. 3.2 and 3.3, 
respectively. We demonstrate how the knowledge of the symmetry group of the problem allows 
one to linearize the complicated nonlinear equations for Lagrangian densities resulting from the 
integrablity conditions. 

3.1 Lagrangian densities of the form / = u x u y g(u t ) 

The integrability conditions (11.2j) imply a single fourth order ODE for g(z), 

g""(g 2 g" - 2g(g') 2 ) - V) 2 + 2^W" + 8(g') 3 g'" - g\g'"f = 0; (3.7) 

to comply with the standard notation, the argument of g is now denoted by z. This equation 
enjoys a remarkable SL(2, R)-invariance inherited from (I2.4p : 

az + 3 „ . „ . 
^=^-y, g = (iz + 5)g; (3.8) 

here a, (3, 7, 6 are arbitrary constants such that aS — d'j = 1. Moreover, there is an obvious 
scaling symmetry g — ► Xg. The equation (13. 7p can be linearized as follows. Introducing h = 
g'/g, we first rewrite it in the form 

h"'(ti - h 2 ) = h 6 - 3h*ti + 9h 2 (h') 2 - 3(h') 3 - Ah 3 h" + (h") 2 ; (3.9) 
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the corresponding symmetry group modifies to 

az + ft 



72 + 5 



h = { riz + 5yh + 1 ( 1 z + 5). (3.10) 



We point out that the same symmetry occurs in the case of the Chazy equation QlJ, p. 342), 
as well as its analogue discussed recently in [2]. The presence of the SL(2, M)-symmetry of this 
type implies the linearizability of the equation under study. One can formulate the following 
general statement which is, in fact, contained in [5]: 

Theorem 3. Any third order ODE of the form F(z, h, h', h", h!") = 0, which is invariant 
under the action of SL(2,M) as specified by A3.10\) , can be linearized by a substitution 

z = —, h = -^-\nw 2 (3-11) 

w 2 dz 

where W\(t) and Wz(t) are two linearly independent solutions of a linear equation d 2 w/dt 2 = 
V(t)w with the Wronskian W normalized asW = W2dw\j dt — w\dwij dt = 1 (the potential V(t) 
depends on the given third order ODE, and can be effectively reconstructed) . 

In particular, the general solution of the equation Ii3. 9\) is given by parametric formulae 
(EX!]) where wi(t) and W2{t) are two linearly independent solutions of the linear equation 
d 2 w/dt 2 = §(cosh~ 2 i)u; with W = 1. 

Proof: 

To establish the first part of the theorem we essentially reproduce the calculation from Sect. 5 
in [5]. Let us consider a linear ODE d 2 w/dt 2 = V(t)w, take two linearly independent solutions 
wi(t), W2(t) with W = 1, and introduce the new dependent and independent variables h,z by 
parametric relations 

Wi d 
z = — , h = — In w 2 - 

w 2 dz 

Using the readily verifiable formulae dt/dz = w\ and h = W2dui2/dt, one obtains the identities 

h' -h 2 = w\ V, 
h" - 6hti + Ah 3 = w & 2 dV/dt 

and 

h'" - 12/i/i" - 6(h') 2 + 48h 2 ti - 24/i 4 = iyf d 2 V/dt 2 - 
here prime denotes differentiation with respect to z. Thus, one arrives at the relations 

T _ (/i"-6/ife / +4fe 3 ) 2 _ (dV/dt) 2 
1 1 — (h'-h 2 ) 3 — V 3 ' 

j _ h'" ~12hh" ~6(h') 2 +A8h 2 h' -24fe 4 _ d 2 V/dt 2 
J 2 — (h'-h 2 ) 2 — V 2 ■ 
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We point out that l\ and I 2 are the simplest second- and third-order differential invariants of 
the action (13.101) whose infinitesimal generators, prolonged to the third jets z, h, h', h", h'", are 
of the form 

Xx = d z , X 2 = zd z - hd h - 2h'd h , - 3h"d h » - Ati"d h ,„, 

X 3 = z 2 d z - (2zh + l)d h - (2h + 4zti)d h , - (Qh' + Qzh")d h „ - (12h" + &zh!")d h ,n] 

notice the standard commutation relations [Xi, X 2 ] = Xi, [Xi,X 3 ] = 2X 2 , [X 2 , X 3 ] = X 3 . One 
can verify that the Lie derivatives of Ii,I 2 with respect to Xi, X 2 , X 3 are indeed zero. Thus, 
any third order ODE which is invariant under the SL(2, R)-action f 1 3 . X 1) . can be represented 
in the form I 2 = F(Ii) where F is an arbitrary function of one variable. The corresponding 

potential V(t) has to satisfy the relation d V yf l = F ^ ^ dV yf^ 

This simple scheme produces some of the well-known equations, for instance, the relation 
I 2 = —24 implies the Chazy equation for h, that is, h'" — 12hh" +18(h') = 0. The corresponding 
potential satisfies the equation d 2 V/dt 2 = —24V 2 . 

Similarly, the choice J 2 = h - 8 results in the ODE h'" = Ahh" - 2(h') 2 + {h " h , 2 _f? 2 which, 
under the substitution h = y/2, coincides with the equation (4.7) from |2J. The potential V 
satisfies the equation Vd 2 V/dt 2 = (dV/dt) 2 - 8V 3 . 

Finally, the relation 7 2 = Ii — 9 coincides with (13. 9p . The corresponding potential V satisfies 
the equation Vd 2 V/dt 2 = (dV/dt) 2 — 9V 3 . It remains to point out that, up to elementary 
equivalence transformations, the general solution of the last equation for V is given by V = 
I cosh -2 1. 

Since the equation d 2 w/dt 2 = |cosh _2 t w is related to the hypergeometric equation s(l — 
s)w ss + (1 — 2s)w s — |w = 0, corresponding to the parameter values a = 1/3, b = 2/3, c = 1, 
by a change of variables s/ (1 — s) = e 2t , we can reformulate the above Theorem as follows: 

Proposition 1. The general solution of the equation A3. 9\) is given by parametric formulae 
A3.11\) where wi(s) and w 2 (s) are two linearly independent solutions of the hypergeometric 



— ->-)H\ - 

Widw 2 /ds = l/(2s(l — s)) 



equation s(l — s)w ss + (1 — 2s)w s — §ii> = with the Wronskian normalized as w 2 dwi/ds — 



As h = g'/g, this immediately implies the following formula for the general solution of (13.71) : 



Proposition 2. The general solution of the equation \3. 7| j is given by parametric formulae 

Wi 

z = — , g = w 2 , 
w 2 

where Wi(s) and w 2 (s) are two linearly independent solutions to the hypergeometric equation 
s(l — s)w ss + (1 — 2s)w s — |w = with the Wronskian normalized as w 2 dwi/ds — Widw 2 /ds = 
l/(2a(l-a)). 
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One can construct the following explicit solution of the equation (13.71) . 

g ( z )= q {a2 - anp2) = l + 6g + 6g 3 + 6g 4 + 12g 7 + ....; (3.12) 

(a,/3)GZ 2 

here q = e 2mz . To get a real- valued solution, one has to restrict z to the imaginary axis. This 
function is known as the Eisenstein series E\$(z). Equivalently, it can be defined by the formula 

g(z) = E h3 (z) = 1 + 6^ Ij^ Xs(d) I q n 

n=l y d\n J 

where %3 denotes the Legendre symbol mod 3 (that is, Xz(.d) — if d = mod 3, Xs{d) = 1 
if d = 1 mod 3, and Xz{d) = — 1 if d = 2 mod 3). The Eisenstein series transforms as 
9( °y Z z +s ) = Xs^Xt 2 + S)g{. z ) under the Hecke congruence subgroup T (3) defined as 

(" ?)er (3)cSL(2,z) if (° f)-(o mod3 - 

It follows that (7(2) is a modular form of weight one and level 3, namely, g(^^§) = {jz + 5)g{z) 
where 

(° ^) eS L(2.Z) and (« *) S (J { ) m „d 3. 
The function g(z) can also be written in the form involving summation over N only, 

/ ^3fc-l ^3fc-2 \ 

^(z) = l-6^( x _ - x _ j. 



Theorem 4. The function g(z) is a solution of the differential equation (3.1). 

Proof: 

Recall that, given a modular form g(z) of weight k, its Rankin-Cohen brackets [g, g}2 and [g, g}^ 
are defined as follows: 

[g,g] 2 = (k + l) (kgg"-(k + l)(g f ) 2 ), 

and 

r 1 (h u_ oVt -l. ^ ^ + 1) - (fc + l)(fc + 3) , ,„ (fe + 2)(fc + 3) 2 
[SSS^ = (fc + 2)(fc + 3) I — — — ^ gg + (g ) 

these are known to be modular forms of weights 2k + 4 and 2A; + 8, respectively (we use the 
normalization of [20J). In our case k = 1, so that we get 

G = [g,g} 2 = 2(gg" - 2(</) 2 ), [g,gU = 2(gg"" - lGg'g'" + 18G/T), 
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weights 6 and 10, respectively. One can verify that, up to a constant multiple, the left hand 
side of the equation (I3.7P can be represented in the form 



7[9,9]i[9,9]2 + [G,G] 2 , 

which shows that it is a modular form (in fact, a cusp form) of weight 16 with respect to the 
same group. To show that this form vanishes identically we recall that the dimension of the 
space of cusp forms of weight 16 and level 3 equals 4, and the order of zero cannot exceed 5. 
Thus, it is sufficient to verify the vanishing of the first five coefficients in the decomposition of 
this form as a power series in q = exp(2iriz). This can be done by a direct calculation. 

Remark. The relation between the modular form (I3.12p and the hypergeometric equation 
from the Proposition 2 can be summarized as follows. Choosing a basis of solutions of the 
hypergeometric equation in the form 

2 10 2 5 57 2 

Wo = 1 H — s H s + wi = W2ins-\ — sH s + 

9 81 ' 9 162 

one obtains parametric equations 

wi , 5 37 2 „ 2 10 2 

z = — =msH — s H s + q = w 2 = 1 H — s H s + .... 

w 2 9 162 ' y 2 9 81 

Solving the first equation for s in the form s = e z + ae 2z + be 3z + ... one arrives at s = 
e z — |e 2z + || e 3z + .... The substitution into the second equation implies the expression for g(z) 
in the form g(z) = 1 + pe z + qe 3z + ... which, up to an appropriate affine transformation of z, 
coincides with ( 13.121) . 



3.2 Lagrangian densities of the form / = g(u x , u y )u t 

As follows from |8], the integrability conditions (11.21) result in a system of five equations ex- 
pressing all fourth order partial derivatives of g(a, b) in terms of its lower order derivatives. In 
symbolic form, this system can be represented as follows: 

d A g = d 3 g^ + 6^-det(dm) + 3^f^det(dn). (3.13) 
h h h 

Here d s g are symmetric differentials of g, the matrices m and n are defined as 



g a 9b 

HI | g a 9aa 9ab | ■ " 

9b 9ab 9bb 



9aa 9ab 
9ab 9bb 



and 

h = -det{m) = g\g aa - 2g a g b g ab + g\g hh . 

The non-degeneracy of the Lagrangian density f(a,b,c) = g(a,b)c is equivalent to the con- 
dition h 7^ 0. One can show that the over-determined system (13. 13j) is in involution, and 
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its solution space is 10-dimensional (indeed, the values of partial derivatives of g up to or- 
der 3 at a point (ao,&o) amount to 10 arbitrary constants). The system (I3.13P is invariant 
under a 10-dimensional group of Lie-point symmetries which consists of arbitrary projective 
transformations of a and b, isomorphic to SL(3, M), along with transformations of the form 
g — > ag + /3, a,/3 = const. The corresponding infinitesimal generators include 

d d 
2 translations : — , — ; 

oa do 

d d d d 

4 linear transformations : a—, b-—, a—, b- 



(3.14) 



9a' da <9o' <%' 

d d d d 

2 projective transformations : a 2 — — h ab— , ab— — h b 2 — ; 

da ab da db 

2 affine transformations of g : 9tt- 

og ag 

We will demonstrate that the existence of this symmetry group allows one to linearize the 
integrability conditions ( 13.131) . For this purpose we consider a linear system of the form 

z xx = Mz x — Izy + Az, 

z xy = -Nz x -Mz y + Bz, (3.15) 

Zyy —— J Z x ~~\~ TV Zy ~\~ C Z , 

where the coefficients J, J, M, N,A,B,C are certain functions of x, y which have to satisfy the 
compatibility conditions resulting from the requirement of consistency of the equations (13.151) : 

A = 2(M 2 + IN) + I y - M x , B = M y + N x + I J - MN, C = 2(N 2 + JM) + J x - N y , 

along with two extra relations involving J, J, M, N only: 

J xx = 2N xy + M yy - 3{N 2 + JM) X - 3NM y + 2JI y + IJ y , 
I yy = 2M xy + N xx - 3(ikf 2 + I N) v - 3MN X + 21 1 x + JI X . 

These relations imply that the space of solutions of the linear system (13.151) is three-dimensional. 
Notice that any involutive second order linear system of the form = T^jZk + QijZ with 
two independent variables x, y can be reduced to the form (I3.15P by a gauge transformation 
z — > <p(x,y)z. This is a standard normalization in the theory of multi-dimensional Schwarzian 
derivatives [17J. It implies that the Wronskian W of any three linearly independent solutions, 

Z\ z 2 Z3 

W = det (zi) x {z 2 ) x (z 3 ) x 

(Zl)y {Z 2 )y (Z 3 ) y 

is constant: W x = W y = 0. Let us choose three linearly independent solutions Zi(x,y) with the 
Wronskian normalized as W = 1, and introduce the new independent variables 

a = — , b — — . 

z 3 z 3 
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Let us consider x and y as functions of a, b: 



x = x(a,b), y = y(a,b). 
A direct calculation shows that these functions satisfy the following nonlinear system: 

X a X bb 2x a X b X ab -|- X^X aa (XaVb Ua-^b) J> 
VlVbb - ^VaVbVab + vtVaa = {^aVb ~ VcFb) 2 I, 

(3.16) 

y a x aa - XaVaa = 3Nx a y 2 a + Jy\ - Ix\ - 3Mx 2 a y a , 

XbUbb - ybXbb = 3Mx\y b + Ix\ - Jyl - 3Nx b y 2 ; 

here /, J, M, N are the same functions of x, y as in (13.151) . The system (13.161) is in involution, 
and its solution space, which is 8-dimensional, possesses a transitive action of SL(3, M): indeed, 
there is an SL (3, IR)-freedom in the choice of a basis 21,22,23. Conversely, one can show that 
any involutive system of four second order PDEs for two functions x(a, b) and y(a, b) which 
is invariant under a transitive projective action of SX(3,R), has the form (I3.16p . and comes 
from a linear system (13.151) . This immediately follows from the equivalent representation of the 
system (I3.16p . 

X a Xbi) 2,X a XijX a b -\- XiX aa , » 

- = J(x,y), 



(XaVb ~ VaXb) 2 

vlvbb - 2y a yby a b + y 2 h y a 



i(x,y), 



(XaVb - VaXb) 2 

(3.17) 

yjxgg + y 2 a x bb + 2x b y b y aa + 2x a y a y bb - 2y a y b x ab - 2(x a y b + x b y a )y ab _ , 

- -- — 6M{x,yj, 

[XaVb ~ VaXbY 

x lVbb + x\y aa + 2x a y a x bb + 2x b y b x aa - 2x a x b y ab - 2(x a y b + x b y a )x ab 



{XaVb ~ VaXb) 2 



-3N(x,y). 



The variables x, y and the left hand sides of (13.171) form a basis of differential invariants of 
the SL(3, IR)-action extended to second order jet space with coordinates a, b, x, y, x a , x b , y a , y b , 
x a a,x ab ,x bb ,y aa ,y ab ,y bb . Thus, any system with the required symmetry properties can be ob- 
tained by expressing the four second order differential invariants as functions of x and y. The 
expressions in the left hand sides of (13.171) are related to the two-dimensional Schwarzian deriva- 
tives rjTj. 

Let us return to the integrability conditions (13.131) . Extending the action of the symme- 
try generators (13.141) to the third order jet space with coordinates a,b, g, g a , g b , g aa , g a b, gbb, 
g aaa , g aa b, g a bb, gbbb according to the prolongation formulae (12.61) . one obtains 10 vector fields on 
a 12-dimensional space; thus, there exist two differential invariants, which we will denote by 
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x and y, respectively (the explicit formulae for x and y in terms of g and its derivatives are 
provided below). We claim that x and y, viewed as functions of a and b, satisfy a system of 
the form (13.161) (equivalently, (I3.17P ). Indeed, the action of the SL(3, R) part of the symmetry 
group on the variables a, b, x, y is exactly the same as for the system (I3.16p . The passage from 
a, b, g to a, b, x, y can be viewed as a factorization of the integrability equations (I3.13P by the 
two-dimensional affine group corresponding to the last two generators (13.141) . which respects 
the action of SL(3, M.). To write down the expressions for the differential invariants x and y we 
introduce the following notation: 

9aaa _ ^_gla 9bbb _ 3 gj h 

yl _ 9a 2 gl 2 4 2 _ (fe 2 g| 4 2 

( 9 O i 2 \9 9 a Ub-> ^ ~ I 9 r> I 9 \ 9 9 a !Jb I 

to&*> - 2 9a9bg a b + g h gaa) 2 [gigbb - 2g a g b g a b + g h 9aaY 



gaab n g aa g a b , Ifi'fe 2 n 9bb9ab . 1 5a 



£/a ^3^3 ^2 __ 9b gb z 5b 

(glgbb - 2g a g b g a b + glg aa ) 2 a b ' (glgbb - 2g a g b g a b + g^g, 1,2 



V - („9 . o . . . .2.. D 5,,///,- * ■,.;». .2 . r> -'A, '//, • 

Moreover, let 

S = - 3V 1 - (Z 2 - 3V 2 ), 5 = Z 1 + Z 2 - V 1 - V 2 . 
Then the invariants x and y can be chosen in the following form: 

[S 3 + 9£S + rS^ 1 - Z 2 )} 2 T 2 - S 2 + 4£(V X - V 2 ) + 4(^ + V 2 ) - 1 

X ~ [£ 2 + 6S + 3] 3 ' y ~ [£ 2 + 65 + 3] 2 ' 

notice that the expressions for x and y are manifestly symmetric under the interchange of 
indices 1 <-> 2. The functions x(a, b) and y(a, b) satisfy a system of the form (I3.16P which 
is invariant under the action of SL(3, R) as specified above. The explicit formulae for the 
coefficients /, J, M, N, as well as the properties of the corresponding linear system (I3.15p . will 
be discussed elsewhere. 



3.3 General case: Lagrangian densities of the form f(u x ,u y ,ut) 

Let us begin by introducing special functions which will appear in the general formula for /. 
The first one is a theta-function of order one [12J (with modular parameter r — e), 

6{z) = ]T(-l) fc exp Uri{kz + ^^e)) , (3-18) 

which is known to satisfy the relations 

6(z + 1) = 9{z), 6(z + e) = - exp(-2mz)B(z)] 

here e = |(1 + Voz). Note that e is a primitive 6th root of unity. In particular, e 3 = — 1 and 
e 2 = e — 1. It is known that 0(0) = 0, and this is the only zero of the function 9{z) modulo 1, e. 
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Moreover, 8'(0) ^ 0. Next we define a function 8(z) which differs from 8(z) by an exponential 
factor: 

~ 6{z) = W) exp (~ 27ri( l ~ \ )z2 ~ ™) 6{z) ' (3,19) 

It can be readily verified that 

0{z + 1) = exp(-^f((2 £ -l) z + e- 2))8(z), 



(3.20) 



8(z + e) = exp(-^((e + l)z + e- 2)))8{z). 

These relations imply 

9{z + a + (3s) = exp f-2ni((^±a + ^/3)z + + £ -^a(5 + ^(3 2 ) ) (:>(: ) 

V o o o o o 



(3.21) 



where a, /5 £ Z. One can also show that 

£(e2) = e9(z), 

which implies 



§(z) = j2 a i z6j+1 -> ( 3 - 22 ) 

j>0 



note that a = 0'(0) = 1. 

Finally, let us define the 'master-density' / by the formula 



(fc,/)ez 2 \o 



which we claim to be a 'generic ' solution of (fOf l For convenience, the arguments of / are 
now denoted by x, y, z. From (13.211) one can derive the following modular properties of /: 

f(x + l,y,z+ (2e - l)x + e + l) = f{x, y, z) + y, 

f{x + e, y, z + {£ + l)x + e + 1) = f{x, y, z) + ey, (3.24) 
f(x, y + 1, z + (2e - l)y + e + 1) = f(x, y, z) + x, 
f(x, y + e, z + (e + l)y + £ + 1) = /(a;, y, 2) + ex. 
Substituting (13.221) into f!3.23j) we obtain an alternative representation for /, 

f(x,y,z)=xy+ a m a n x 6m+ V n+1 ^-d) 6 ( m+n )exp^(A; 2 -^ + / 2 )^ = 

(fc,/)eZ 2 \0, m,n>0 ^ ^ ' ' 



1 We did not prove that / satisfies the system (11.2[) . although computer calculations support this claim. 
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m,n>0 

where 

9n(z) = - d ) 6n + ( £k ~ Z ) 6n ) 6X P ~ H + ■ 

It is known [13] that g n (z) is a modular form, namely, g n ( ^^§ ) — (7-2 + ^) 6n+1 fi , n(- 2 ) where 
f ^ ^ j G SX(2, Z) and ^ ^ ^ ^ = ^ ^ ^ ^mod 3. Indeed, coefficients at the exponents 

are harmonic polynomials with respect to the quadratic form k 2 — kl + I 2 . This gives us the 
following modular property of f(x,y, z) with respect to the same group: 

t ( x y az + (3\ 1 

/ ?) f> F = J\ x ^y^ z ) F- (3.25) 

\^z + o r )z + o r yz + oJ 72 + 

Note that g(z) = go(z) is a solution of the differential equation (13. 7p and is given by (I3.12p 
(up to a rescaling of 2). Functions g n (<2) can be represented as rational differential functions in 
g(z), for example, g\{z) = g(z) 2 g"(z) — 2g(z)g'(z) 2 = ^g[g,g} 2 where [g,g}2 is the Rankin-Cohen 
bracket (see the proof of Theorem 4). Transformations (13.241) and (I3.25|) generate a discrete 
subgroup r C SX(4, C) which plays the role of the modular group for f(x, y, z). The subgroup 
r has the following algebraic description. Let D C C 3 be a domain in C 3 defined by 

2 

D = {(x,y,z) G C 3 ; \x\ 2 + \y\ 2 < -pQz}. 

V3 

Note that D is a complex hyperbolic ball. One can check that the series (I3.23P converges exactly 
in the domain D. Let G C SX(4, C) be a group defined by 

G = {A G SL{4, C); A J A* = J} 

where 



J 



( 1 \ 
10 
-y/3i 
\ V3i J 

and A* stands for the Hermitian conjugate of A. One can check that the complex hyperbolic 
ball D is an orbit of G under its standard projective action on CP 3 : if A = (a^) G G and 
(x, y, z) — (x : y : z : 1) G CP 3 , then 



A(x,y,z) 



a u x + a 12 y + a 13 z + a u a 21 x + a 22 y + a 23 z + a 24 a 31 x + a 32 y + a 33 z + a 34 
dux + a^ 2 y + a i3 z + a 44 ' a 44 x + a i2 y + a^ 3 z + a 44 ' a 41 x + a^ 2 y + a 4 3Z + a 44 



In our context, the domain D and the group G play a role similar to that of the upper half plane 
and its automorphism group SL(2,M) in the classical theory of modular forms. Let fcGbe 
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A = 



mod (1 + e). 



a discrete subgroup of G consisting of matrices A = (ay) G G with the following properties: 
ctij G Z[e] and 

/ 1 a 13 au \ 
1 a 23 a 2 A 
1 
\ a 43 1 / 

Here Z[e] = {m + en; n,m G Z}, and a = b mod(l + e) for a, b G Z[e] means |=| G Z[e]. 
Note that all matrices corresponding to the transformations (13.241) and (13.251) belong to T. We 
conjecture that the group T is generated by these transformations. 

Remark 1: limiting cases. Lagrangian densities of the form xg(y,z) can be obtained as 
linit^o ^ tx ^ ,z ^ from f(x, y, z) as defined by (I3.23p . This gives the function g(y, z) in the form 

g(v,z)=y+ J2 m k Z fi V) ex P (ifV ~ ** + = E a ny 6n+1 9n(z). (3.26) 

This function is defined on the domain {(y, z) G C 2 ; |y| 2 < -^^ z } and satisfies the following 
modular properties: 

g(y + 1, z + (2e - l)y + e + 1) = 2) + 1, 

<?(y + e, z + (e + l)y + e + 1) = g(y, z) + e, (3.27) 

— rjf > — rr = 9(v> z >- 

\72 + 7^ + / 

Similarly, Lagrangian densities of the form xyg(z) can be obtained as lim^o ^ tx 'l y ' z ^ . This 
brings us back to the modular form g(z) discussed in Sect. 3.1. 

Remark 2. Computer experiments show that solutions of the system f 1 1 . 2 1) can also be 
sought in the form of a power series, 

f{x,y,z)= c ijk x 6i +Y j+1 z 6k+1 . 
i,j,k>a 

Moreover, = ai(ija k b i+ j +k where are the same as in (I3.22p . and 6j is yet another sequence 
of complex numbers. Taking a limit we obtain 

g(y,z)= J2^ k b 1+ ky 6j+1 z 6h+1 

for densities of the form f(x, y, z) = xg(y, z), and 

g(z) = J2^b k z 6k+1 

k>0 

for densities of the form f(x, y, z) = xyg(z). We point out that integrable densities of this type 
differ from the solutions constructed above by appropriate transformations from the equivalence 
group. 



15 



Remark 3. The action of the 20-dimensional equivalence group on the 23-dimensional 
space of third order jets of the function f(a,b,c) possesses three differential invariants which 
we denote by x, y, z. Viewed as functions of a, b, c, these invariants satisfy a nonlinear system 
which possesses a transitive action of SL(4, R). This system can be linearised following the 
procedure outlined in Sect. 3.2. We plan to report the details elsewhere. 

Remark 4. There exist a number of examples of integrable Lagrangian densities expressible 
in terms of elementary functions. One can mention, e.g., the four polynomial Lagrangians 
classified in [7]: 

f=U X UyU t , f = U 2 x Uy + UyU t , f = /3 + ~ U X U t 

and 

f = u 4 x + 2u 2 x u t - u x u y - u\. 

It would be interesting to explicitly demonstrate how these (and other) examples can be ob- 
tained as degenerations of the 'master-Lagrangian' constructed in Section 3.2, and to describe 
singular orbits of lower dimensions. 
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